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The Gardiner's phonon presented for a particle-number conserving approximation method to de- 
scribe the dynamics of Bose-Einstein Condesation (BEC) (CW. Gardiner, Phys. Rev. A 56, 1414 
(1997)) is shown to be a physical realization of the g-deformed boson, which was abstractly devel- 
oped in quantum group theory. On this observation, the coherent output of BEC atoms driven by a 
radio frequency (r.f) field is analyzed in the viewpoint of a g-deformed Fock space. It is illustrated 
that the q-deformation of bosonic commutation relation corresponds to the non-ideal BEC with the 
finite particle number N of condensated atoms. Up to order 1/N, the coherent output state of the 
untrapped atoms minimizes the uncertainty relation like a coherent state does in the ideal case of 
BEC that N approaches infinity or q = 1 — 2/N approaches one. 



I. INTRODUCTION 

Since the introduction of the concept of q-deformed bo- 
son by different authors independently in the develop- 
ment of the quantum group theory , there have been 
many efforts of finding its physical realizations |^||, e.g., 
fittings the deformed spectra of rotation and oscillation 
for molecules and nuclei In our opinion, those in- 

vestigations can be regarded as merely phenomenological 
because a g-deformed structure is postulated in advance 
without giving it a microscopic mechanism. In this paper 
it will be shown that a physical and natural realization 
of the g-deformed boson is provided by the phonon oper- 
ators, which was presented recently by Gardiner JhJ for 
the description of Bose-Einstein condensation (BEC). 

As well known, to deal with the dynamics of Bose- 
Einstein condensated gas the Bogoliubov approximation 
in quantum many-body theory is an efficient approach, 
in which the creation and annihilation operators for con- 
densated atoms are substituted by c-numbers. One short- 
coming of this method is that the total atomic particle- 
number may not be conserved after the approximation. 
Or a symmetry may be broken. To remedy this default, 
Gardiner suggested a modified Bogoliubov approxima- 
tion by introducing phonon operators which conserve the 
total atomic particle number N and obey the bosonic 
commutation relations in case of N—f oo. In this sense 
this phonon operator approach gives an elegant infinite 
atomic particle number approximation theory for BEC 
taking into account the conservation of the total atomic 
number. 

What will be investigated here is the case that the 
total atomic particle number TV is very large but not 
infinite. That is, we shall consider the effects of order 
o(l/N). And we shall focus on an algebraic method of 
treating the effects of finite condensated particle number 
in atomic BEC. As it turns out, the commutation rela- 



tions for the Gardiner's phonon operators will no longer 
obey the commutation relation of the Heisenberg-Weyl 
algebra but the g-deformed bosonic commutation rela- 
tion 

[Ml, = btf-qtfb = l, (1) 

where the deformation constant q depends on the total 
atomic particle number. By using this g-deformed bo- 
son algebra we will discuss the effects of the finiteness of 
the particle number on the output of MIT experiment of 
BEC. 

This paper is organized as follows. In sec. II after 
the reformulation of the Gardiner's modified Bogoliubov 
method for the two-level atoms, the g-deformed commu- 
tation relation for the Gardiner's phonon operator is de- 
rived in case of large but finite particle number. The 
excitation of the Gardiner's phonons is also considered 
for small N. In sec. Ill, the Hamiltonian describing the 
coherent output in MIT experiment for atomic laser Jp| 
is discussed in terms of the g-deformed algebra. When 
q = 1 or N — ► oo it describes the coherent output of 
BEC from a vacuum to a coherent states. For general q, 
or finite but large condensated atomic particle number, 
it rusults in a nonlinear dynamic Hciscnbcrg equation 
for the g-deformed bosonic operators. Its solution at the 
first order of 1/N is analyzed. In sec. IV dressed phonon, 
a generalization of the Gardiner's phonon in case of quan- 
tized radiation field, is disscussed and it satisfies also a 
g-deformed bosonic algebra for a large phonton number 
and atomic particle number. Finally there are some con- 
cluding remarks in the last section. 

II. g-DEFORMED BOSONIC ALGEBRA FOR 
THE GARDINER'S PHONON 

Gardiner's starting point to introduce the phonon op- 
erators is to consider a system of the weakly interacting 
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Bose gas. Without losing generality, we consider for the 
moment a classical radiation field interacting with two- 
level atoms with b\ and 6 e denoting the creation and an- 
nihilation operators for the atoms in the untrapped state 
and &J and b g for the creation and annihilation operators 
of the atoms in the trapped ground state. The simplified 
Hamiltonian for the MIT experiment reads 



H = hu e b\b e + fig(t)b\b g + h.c. 



(2) 



where g{t) is a time-dependent coupling coefficient for the 
classical sweeping r.f. field coupled to those two states 
with level difference fiu> e - Note that the total atomic 
particle number N — b\b e + b' g b g is conserved. For con- 
venience we define 77 = 1/iV for large particle number. 

In the thermodynamical limit N — > oo, the Bogoli- 
ubov approximation is usually applied, in which the lad- 
der operators b g ,b g of the ground state are replaced by 
a c-number >/N c , where N c is the average number of 
the ininital condensated atoms. As a result Hamiltonian 
Eq. (||) becomes a forced harmonic oscillator with the ex- 
ternal force term Ti\fW~ c g{t)b\ + h.c. Obviously, it will 
drive the atom in the untrapped excited state from a 
vacuum state to a coherent state and thus lead to a co- 
herent output of BEC atoms in the propagating mode . 
In this case the condensated atoms are initially described 
by a coherent state JT3| . 

However, this apporoximation destroyes a symmetry 
of the Hamiltonian Eq.(Q), i.e., the conservation of the 
total particle number is violated because of [N, ^ 0. 
According to Gardiner, the phonon operators are defined 
as: 
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These operators act invariantly on the subspace V N 
spanned by bases \N; n) = \N — n,n) (n — 0,1, ... , N), 
where Fock sates 



\m, n) 
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\lm\n\ 



bt m b q n \0) 



with m, n = 0, 1, 2, . . . spann the Fock space Hib of a two 
mode boson. 

A straightforward calculation leads to the following 
commutation relation between the phonon operaotr and 
it Hermitian conjugate: 



M] = l--btb e = f(b^b; V ), 



(4) 



with f(x; rf) = yl + 2(1 — 2x)i] + rf — 77. Keeping only 
the lowest order of 77 for a very large total atomic particle 
number, the commutator above becomes 



[b, 6 f ] = 1 - 2776*6 



(5a) 



(5b) 



with q = 1 — 2r). This is exactly a typical g-deformed 
commutation relation. As N — ► 00 or q — ► 1, the usual 
commutation relation of Heisenberg-Weyl algebra is re- 
gained. 

Since the Gardiner's phonon excitation for BEC can 
be mathematically understood in terms of q-deformed bo- 
son system, the algebraic method previously developed in 
quantum group theory can be of use for the construction 
of a g-deformed Fock space. The existence of a vacuum 
state |0) such that 6|0) =0 is guaranteed by the funda- 
mental structure of the (/-deformed boson algebra p| . 
Starting with the vacuum, from relation 



66 r = q n b r b+ (ri)6 tj 



(6) 



where the g-deformed c-number (n) = n — n(n — 1)77) 
is up to the first order of 77 and approaches to n when 
N — > 00, the g-deformed Fock space for the Gardiner's 
phonon excitation of BEC is spanned by bases 
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=6 tn |0) 



with (n)\ = (n) (n - 1) ■■•(!) (77 = 0,1, ...,N). The 
action of operators 6 and 6* are then represented on the 
deformed Fock space as: 



6^77) = y/(n + l)\n + I) , b\n) = y/ffi\n - 1) (8) 

the actions of the creation and annihilation operators 6 
and 6* coincide with their direct actions on the subspace 
V N of the Fock space of a two-mode boson on which the 
total particle number is conserved. 

Actually, the g-deformed vacuum state is physically a 
two-mode boson state \N; 0) and it is trivially annihi- 
lated by the phonon operator 6. So the vacuum state 
can be regarded as a case that all bosons condensed in 
the ground state and tends to form a BEC. By map- 
ping back to two-mode boson space, the g-deformed Fock 
state \n) ~ \N;n) implies a phonon excitation that n 
bosons are created in the propagating mode by destroy- 
ing n bosons in trapped mode. This kind of fundamental 
process of phonon excitation essentially guaranteed the 
particle number conservation. 

In the above discussion about the Gardiner's phonon 
excitation, we have linearized commutator h = /(6*6; 77) 
so that a g-deformed commutation rule was obtained. 
Essentially this linearization establishes a physical real- 
ization of the g-deformed algebra. However, if the total 
particle number N is not large enough, then h can not 
be approximated by a linear function. From the commu- 
tation relations between h and 6, 6* 



M] 



[h,b] = -b, 



(9) 



we see that the algebra of Gardiner's operators is a rescal- 
ing of algebra su(2) with factor N. When N — > 00, it 
approaches the Heisenberg-Weyl algebra because h ap- 
proaches identity. Its representation space is a finite di- 
mensional space is found to be spanned by bases given in 
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Eq.(0) with ( n) = n — n(n — l)r) being exact. It is inter- 
esting to notice that (n) coincide with the g-deformed 
number in the first order of rj. The action of those 
creation and annihlation operators have the same form. 
On this finite dimensional space we have b\0) — and 
h\0) = |0) together with tf\N) =0. It is worthy to point 
out that the last equality can not be composed to the 
usual Heisenberg-Weyl algebra because such a constrain 
will make it is impossible for the creation operator to 
be the Hermitian conjugate to the annihilation operator. 
Physically the usual boson algebra and its g-deformation 
with q being not a root of unity describe systems with- 
out particle number conservation such as the black body 
radiation, but this finite dimensional representation pro- 
vide us with a suitable algebraic framework to treat the 
dynamic process where the particle number is conserved. 



III. THE COHERENT OUTPUT OF BEC AS A 
Q-DEFORMED BOSONIC EXCITATION 

In this section we shall apply the concept of the q- 
deformed boson algebra for the Gardiner's phonon exci- 
tation to the theoretical analysis of the MIT experiment 
for the coherent output of the BEC atoms from a trap. 

By denoting N e = b\b e , in terms of the deformed boson 
operators the Hamiltonian can be rewritten as: 



H q = hui e N e 



hVN(g(t)tf +h.c. 



(10) 



For an ideal BEC with infinite atoms condensated in the 
trapped state, the g-deformed bosonic commutation rela- 
tion becomes the usual commutation relation [6, b^\ = 1 
of HW algebra. Together with commutator [N e ,b] = —b, 
a linear Heisenberg equation 



b = —iu> e b — ifJ.{t) 



(11) 



is acquired where [i(t) = g(t)y/N is finite as N tends to 
infinity. Its solution can be easily found to be 



b{t) = 6(0)e" 



f3(t) 



(12) 



with a non-vanishing expectation value in the vacuum 
state: 



m = (o\b(t)\o) 



-lie 
'o 



= (t -*' ) n(t')dt' '. (13) 



Therefore when the BEC is initially at a Fock state, one 
obtains a coherent state output of the untrapped atoms. 
This is just a direct manifestation of the quantum coher- 
ence through a factorization (w (t) b (t)) — (b* (t)) (b(t)) . In 
fact for infinite N, the average phonon number (b*(t)b(t)) 
approach the average atomic number b\b e of outputting 
in the untrapped mode. 

Further, to consider the quantum coherence of the out- 
putting atoms represented by phonon operators b(t), we 



evaluate the variances of the Hermitian quadrature phase 
operators [M 



X 1 (t) = 



V2 • 



X 2 (t) 



iV2 



(14) 



One can easily find that 



A 2 X x , 2 (t) 



\Xl 2 \0) -(0|AX lj2 |0) 2 = - (15) 



are time-independent constants that minimize the uncer- 
tainty product A 2 Xi(t)A 2 X 2 (t) = 1/4. It reflects that 
the outputted atoms are just in a coherent state. This 
judgements for the case with infinite N can be similarly 
used to analyses the influence of g-deformation with finite 
atomic number on the dynamics of outputting BEC. 

In practice condensated atomic number N cannot be 
infinite and therefore the Bogoliubov approximation can 
not work well and so an approximation method with 
higher accuracy that for the generic Bogoliubov ap- 
proximation or its Gardiners modification. And the q- 
deformed boson algebra disscussed in the previous sec- 
tion provides us with a algebric structure of treating this 
correction of finite atomic particle number. 

Up to the approximation of order 77, as discussed previ- 
ously, the dynamics of Gardiner's phonon excitation for 
coherent output of BEC is determined by the g-deformcd 
boson commutation relation. The corresponding Heisen- 
berg equation 

b(t) = -iuj e b(t) - ifi(t) + 2ir ) ^{t)b if {t)b(t) (16) 

is then a nonlinear one. It is should be remarked here 
that that the effective coupling fi(t) = g(t)VN is finite 
in the thermodynamical limit of an infinite number of 
atoms in an infinite volume but with the density fixed, 
i.e., n = limjv y— >oo (N/V) is finite. This is because 
the coupling constant g(t) between the r.f. pulse and 
atoms must be proportional to the inverse of the square 
root of t he th e effective mode volume V. So we have 
fj,(t) oc yjN/V, the square root of the density of the con- 
densated atoms in the thermodynamical limit. For this 
reason the nonlinear terms in Eq.(16) can be handled 
as a perturbation in comparison with other terms. Oth- 
erwise, if /x(t) were not fixed as N approaches infinity, 
the non-liner term should not be regarded as a perturba- 
tion due to its square amplification of perturbation. The 
thermodynamical limit is our key point to deal with the 
dynamics of q-deformation for the coherent output of the 
BEC atoms with a perturbation method. 

Becasue of the nature of the deformed bonon, in gen- 
eral, it is only neccissary to solve this nonlinear Heisen- 
berg equation Eq. ( |l6| ) according to the first order of rj by 
perturbation. In fact if the annihilation operator at time 
t is expanded as 



b(t) = b e- 



(3(t) + vh(t) 



(17) 



with 60 = b(0), one obtains immediately 
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h(t) = £(t)b + ^(ty^bi - W{t)b% + a(t), (18a) 

a(t)=2i[ e-^-^fMit'^Pit'^dt', (18b) 
Jo 

£(t) = 2i[ e-^VWC*')^'- (18c) 
Jo 

As a result the average outputted atomic number 
(N e (t)) in the untrapped excited state with finite con- 
denstated BEC particle number correction is given by 

(N e (t)) = \m\ 2 +v(\mf + ^e(p(t)a(tr)) ■ (19) 

The quantum fluctuation of the quadrature phase op- 
erators Xi(t) and X2(t) in the final state can also be 
explicitly determined as: 

A 2 *i = \-fl [W)\ 2 + Pit) + P* 2 (t)) , (20a) 
A 2 X 2 = 1 - „ (\(3(t)\ 2 ~ P 2 (t) ~ P* 2 {t)) ■ (20b) 

Keeping the first order of 77, one obains immediately the 
follwing relation 

A 2 X 1 A 2 X 2 = \~ r l l/3(*)| 2 • (21) 

The term of 77 reflects the devotion from the thermody- 
namical limit as the effect of the gr-deformation with finite 
condensated particle number. 

Since the commutation relation of quadrature phase 
operators Xi(t) and X2(t) is deformed into 

[X 1 (t),X 2 (t)}=i(l-2r ] b\t)b(t)), (22) 

the corresponding Heisenberg's uncertainty relationship 
of those quadratures in the output taking account of fi- 
nite atomic number should be: 

A 2 X 1 A 2 X 2 > 1 - Tj<0|6t( t )6(t)|0). (23) 

Considering (0\tf (t)b(t)\0) = \f3(t)\ 2 + 0(77) the coherent 
output state of the untrapped atoms minimizes the un- 
certainty relation. Because, as we known, the squeezed 
state and coherent state can be defined to be states 
minimize the uncertainty relation between the position 
and momentum, the output state of the atomic lase can 
be likened to a kind of ^-deformed coherent state or 
even a squeezed state because fluctuations Eq.(26a) and 
Eq.(26b) are inequal. 

IV. DRESSED PHONON OPERATOR 

When the r.f field is a quantized one which is described 
by a single-mode boson whose annihilation operator is 
denoted as a, we consider the following Hamiltonian: 

H d = uo e N e + uj g N g + lj N + g{a!b\b e + ab g bt), (24) 



where No = a^a is the photon number operator and g is a 
real coupling constant. Obviously there is a conservative 
quantity A = No + iV e in addition to N. We now define 
the dressed phonon operators as 

B = -^=a)b\b e , B f = —^ab q b\. (25) 

Vna 9 Vna 9 e K ' 

For convenience we denote 770 = 1/A for large A. The 
dressed phonon operators satisfy the following commuta- 
tion relation: 

[B, St] = 1 - 2(770 + V )N e + 77 77(3iV L 2 - N e ). (26) 

In terms of those dressed phonon operators the Hamilto- 
nian Eq. (pi|) can be rewritten as 

H d =u g N + u]uA + uj s N e + ii d {B + Bt). (27) 

where [i d = gvNA and ujs = uj e — uj g — luq. 

As N (or A) tends to infinity the dressed phonon op- 
erators B, B< and N e will give a representation of al- 
gebra su(2). As both N and A tends infinity one ob- 
tains Heisenberg-Weyl algebra [B,B*] = 1. When we 
keep the first order in 0(77) and 0(770), and considering 
B^B = N e + 0(77 + 7/0), we obtain a g-deformed boson 
algebra 

[B,B\ = 1 (28) 

with the deforming constant given by q d = 1 — 2(77 + 7/0). 

Within this kind of approximation, the Heisenberg 
equation of motion satisfied by the dressed phonon is 
given by: 

B = -ilusB - ipd + 2i( Vo + rj)n d BtB. (29) 

Here the coupling constant g is proportional to the in- 
verse of the square root of volume V which makes /i d 
finite as N and A app roaches infinity. Because the evolu- 
tion equation Eq.([28|) is the same with Heisenberg equa- 
tion Eq. (|l6|) , all the results discussed in the previous sec- 
tion should also be valid in case of quantized r.f. field. 

V. REMARKS 

To conclude this paper, we treat the same problem 
in the Schrodinger picture. Without the effect of q- 
deformation caused by finite N, the dynamic process 
can be described as a factorized evolution of the wave 
function |Tq| . Suppose the initial state is \<j>(0)) — \a = 
yN) g <8> |0) e i.e., all atoms occupied the trapped groud 
state with a coherent state while there is no atom in 
the propagating mode. Driven by a r.f. field g(t) = 
geicp(iujft) (g oc y/l/V) with a fixed frequency Wf, the 
system will have exactly a factorized evolution wave func- 
tion of the form \<j>(t)) = \a) g ® \(3{t)) e even for a finite 
N. Such a factorization structure of wave function first 
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pointed out by Mewes et al is very crucial to realize their 
experiment, and its dynamical origins mainly depend on 
the linearity of the original coupling system of b e and 
b g . The coherent state of the second component |/?(t)) e 
of \<t>(t)) means a coherent output of atoms in propa- 
gating mode in the MIT experiment. The influence of 
finite N on Schroedinger wave funct ion is manifested by 
the outputting amplitude \(3(t) \ cx y/N/V\ sin(fit)| where 
Q = yj(uj - u f ) 2 /4 + g 2 . Because the factor y/N/V 
is finite and O tends to \u — w/| in the thermodynam- 
ics limit, the results from the Bogoliubov approximation 
mentioned in Sec. II is just recovered by this exact solu- 
tion. 

There is another nonlinear effect different from that re- 
sulting from the g-deformation with finite N. In presence 
of interatomic interaction in the trapped or untrapped 
states, which roughly is of the form b\bjbkh, the corre- 
sponding Heiscnbcrg equations are no longer linear and 
thus the factorization structure will be broken down. It 
is noticed that such a non-linearity is essentially different 
from that resulting from the q-deformation with finite N 
in this paper. The later originates from the subsystem 
isolation of a large system. 
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